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Abstract 



The conductance of a normal-metal mesoscopic system in proximity to 

superconducting electrode(s) is calculated. The normal-metal part may have a 

general geometry, and is described as a "circuit" with "leads" and "junctions". 

The junctions are each ascribed a scattering matrix which is averaged over 

the circular orthogonal ensemble, using recently-developed techniques. The 

results for the electrical conductance reproduce and extend Nazarov's circuit 

theory, thus bridging between the scattering and the bulk approaches. The 

method is also applied to the heat conductance. 
PACS numbers: 74.50.+r, 74.80.Fp, 73.23.Ps 
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Mesoscopic normal-metal-superconductor (N-S) proximity effects are often described 
theoretically using the Usadel equations [|TJ - the partial differential equations of inhomo- 
geneous (dirty) superconductivity. However, the ohmic resistance of the normal-metal part 
of the structure || can also be found from the "circuit theory of Andreev conductance" 
- a set of algebraic relations recently derived from the Usadel equations by Yuli Nazarov 
|| (the Usadel equations themselves are derived from the Eilenberger equations of clean 
superconductors, which ultimately derive from the Keldysh technique). Here we pursue an 
alternative approach, which associates a random scattering matrix with each of the junctions 
in the "circuit" . These matrices are to be averaged over one of the circular ensembles of ran- 
dom matrix theory, employing diagrammatic techniques which were developed specifically 
for such purposes 

Electron and hole excitations can enter the normal-metal part of such N-S structures 
(the "circuit") from one of the normal electrodes (N), but at low excitation energies they 
can not enter into the superconducting electrodes (S) because of the gap. Instead, Andreev 
reflections may occur — an electron may evolve into a Cooper pair in S, leaving behind 
a hole in the Fermi sea of the normal-metal part |J, with a wavevector related to that of 
the impinging electron by time-reversal. We will assume for simplicity, as did Nazarov [Q, 
that the applied voltages and the ambient temperatures are very small, so that the energy 
dependence of the scattering properties of the structure can be neglected, and the quantum 
mechanical propagation is fully coherent. We also assume the absence of a magnetic field 
(except for Aharonov-Bohm phases), the absence of electron-electron interactions in the 
normal part of the structure, and unbroken time-reversal symmetry and spin degeneracy. 

Unaveraged quantum mechanical transport in a wide variety of structures may be de- 
scribed |7]] using two unitary matrices: Go represents the "non-local" transport in the 
"leads", and S represents random scattering within each of the "junctions", as well as 
Andreev scattering at the external superconducting electrodes. We denote the number of 
internal junctions in the system by K, and the number of junctions and external electrodes 
by M (we have at least one N and one S electrode, so M > K + 2; see Fig. la). The 
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number of transverse modes in the lead connecting junctions k and m is denoted by Wk, m , 
so that iVfc = J2 m Wk t m is the size of the scattering matrix of the kth junction. The A^s 
and the (nonzero) W4, m s are treated as large parmeters, and the conductances are evaluated 
to leading order — weak localization and conductance fluctuation corrections will not be 
considered here. The matrices G and S are of size 2N tot x 2N tot , where N tot = Y^n=i N m 
is the total number of modes, and the factor of 2 represents the electron-hole subspace. Go 
is non-diagonal in the M x M block structure, but diagonal in the electron-hole space. S 
is block-diagonal, with K diagonal blocks consisting of a random scattering matrix for the 
electron-electron sub-block, and the complex conjugate of that matrix for the hole-hole sub- 
block. The remaining diagonal blocks of S either vanish, for an external normal electrode, or 
are off-diagonal in the electron-hole space: the electron-hole Andreev-scattering sub-block 
for a superconducting electrode m is equal to iexp(— ixm) (times a unit N m x N m matrix), 
and the hole-electron sub-block is iexp(ix m )- Here the x m s are the given, time-independent 
phases of the superconducting order parameter in m. The extra factors of i lead to destruc- 
tive interference after two consecutive Andreev scatterings, and are responsible for, e.g., 
proximity-induced gap-like features in the normal part of the structure. 

We now average the K random-scattering diagonal blocks of S over the appropriate 
matrix ensemble — the circular orthogonal ensemble (COE), justifying our expressions both 
physically and technically (see Ref. j|] for details of the latter). Diagrams for the conductance 
may be drawn as in Fig. lb, with matrix multiplication implied along the bottom and the 
top of the diagram. The matrices along the top line are hermitean conjugated, so that 
the diagrams represent probabilities, or amplitudes multiplied by complex conjugates of 
amplitudes. Averaging is evaluated to leading order by connecting the random elements of S, 
represented as double lines, in all possible ways in the plane of the diagram without crossing 
lines. These connections or couplings (Fig. lc) join together n occurrences of the kth random 
scattering matrix and n occurrences of its hermitean (or complex) conjugate, in alternating 
order; the corresponding weights are (— l) n " 1 c n /N^ n with c n = (2n)!/2(2n — l)(n!) 2 (these 
are deduced in Refs. f§] and || from unitarity). Each 2n-fold coupling divides the plane of 
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the diagram into 2n regions, with matrix multiplication performed independently around the 
periphery of each region. The blocks of S which are not averaged over — those representing 
Andreev scattering at the superconducting electrodes — are drawn as black semicircles. 

The averaged propagator G = (Go + GoSGo + . . .) is written as G = Go/ (I— EGr ), where 
the "one particle irreducible part" , or "self-energy" E represents the "average" scattering by 
the junctions (Fig. Id). The structure of E is particularly simple, because (a) both normal 
scattering and off-diagonal (in mode space) electron-hole scattering average to zero due to 
the random phases involved, and (b) the diagonal electron-hole scattering elements of E must 
be equal to each other within each of the K random-scattering blocks, due to ergodicity in 
each junction. We denote the electron- hole averaged Andreev scattering amplitude by iff., 
where is a complex number in the unit circle. The hole-electron amplitude is then if£. 
Direct Andreev scattering (black semicircles) is also conveniently included in E: for the 
superconducting electrodes f m = exp(— ixm), an d for the normal ones f m = 0. 

The behavior of the amplitudes is similar to that of voltages — if a certain junction 
is connected by a lead with a high conductance to an electrode with a given value of f m , 
then the value of fj. in that junction will tend towards f m . These amplitudes are physically 
observable in the following sense: if we were to inject electrons directly into the fcth junction 
from an additional normal electrode which is relatively weakly coupled to that junction - 
a "noninvasive voltage probe" — we would find a strong beam of holes being retro-reflected 
into that voltage probe, with an intensity of \fk\ 2 times that of the electron beam. The fact 
that these amplitudes are of order 1, rather than of order 1/Nj., has been referred to in the 
literature as the "giant (Andreev) backscattering peak" ||. It can be understood on the 
basis of the analogy between Andreev reflections and phase- conjugating mirrors in optics. 

For the specific case of ideal leads, G = G + GoSGo + ... is a simple sum over multiple 
Andreev reflections. For example, G^ h&ke j = if m +if m if%if m + . . ■ = if m /0- + fmfk), where 
j is the index of one of the Wk, m modes in the lead connecting k and m (the superscript 
m has been added to emphasize this; kh and ke denote the kth diagonal electron-hole sub- 
block). On the other hand, the diagrams for E express fk as a sum of powers of and a* k) 

4 



where idk = (1/Nk) Ylf=i Gkh,j-,ke,j is a trace of a sub-block of G. When the explicit values 
of the used and the relationship between au and fk is inverted, it reads simply 

This is consistent with the physical requirement that fk = f m if all the / m s in the sum are 
equal to each other. Eq. ([!]) may be rewritten, using J2 m ^k,m = N k , as an expression of 
"spectral current conservation": 

5>*.m = ° ! Jim = 2 ^ {^+\M 2 )GT hM ~ fk) = W k , m 2 J ™ Jk , (2) 

rn j=l ' JmJk 

where Jk, m denotes the dimensionless (complex) "spectral current" in the (k, m) lead. 

The "two-particle irreducible vertex" , V is represented diagrammatically in Fig. le. It is 
also block- diagonal, with all of the elements of the kih. block given by — — ( 1 ' 



N k {l-\f k \ 2 ) V-l/fcl 2 1 

(each of the four sub-blocks are full matrices with identical elements). This expression is the 
only one consistent with unitarity or current conservation . The two-particle propagator, 
or "diffuson" is given by D = T + TtT + . . . , where the elements of t are equal to the absolute 
squares of the corresponding elements of G. Rather than evaluating the matrix D directly, 
as would be necessary in order to express the conductances in the usual scattering approach 
|?|,|T0|j, we define occupation probabilities v for the electron and hole modes leaving any of 



the junctions or electrodes. For example, tW(/i)j(e) = 1/ (1 + exp[(e ± eV m ) /kBT m ]j for an 
external normal electrode m with an applied voltage V m and a temperature T m (here e is 
the excitation energy relative to the Fermi surface of the superconducting electrodes). The 
values of v for modes leaving internal junctions of the circuit are given by v — Ttv, where 
the probabilities v on the right hand side represent the neighbouring junctions or electrodes, 
and the factors of t and T represent transport into and inside the junction, respectively. It 
is convenient to rewrite this as 

D~ l v = (T- 1 -t)v = 0, (3) 
with the equality required only for the K internal junctions. 



As we assume that no inelastic processes occur in the structure, the electric current / 
and the quasiparticle current Q are conserved differentially at each junction k: 

EW e ) = ° ; E<2*,m(e) = o. (4) 

m m 

The total electric current, Ik, m = / °° de I^ m {e), and heat current, Qk, m = Jo°° e ^ e Qk,m( e ), 
are consequently conserved as well. For ideal leads, Eq. (|3|) implies that || 
T 2e 2 (l + |/ m H(l + l/ fc P) 

& |l + /m/fc| 2 

Q fc , m = ^ fc ,m (1 "l{" |2 l (1 ~ l/fc|2) (^ " ^n) , (6) 



(1- 


fm 
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1 + fmfk 
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where the voltages V m and the excitation-energy densities U m are defined by ]TT 



TA 1 Z" 00 , ^/im ^erra rr Z" 00 7 ^/im "era /_\ 

V m = - de ; U m = dee — — . 7 

eJo 1 + l/ml 2 i-l/ml 

These definitions are reasonable because |/ m | = at the N electrodes. At the S electrodes, 
both the voltages V m and the heat currents Qk, m vanish. The net heat current into an 
external normal electrode m is Qk, m — Ik^Vm and not just Qk, m , because of the shift in the 
Fermi level. This completes the derivation for ideal leads — the conservation of the currents 
of Eqs. (H), (|5|) and @ at each junction determines the Andreev amplitudes fk, the voltages 
Vk and the energy densities Uk self-consistently. Although seemingly local, this scheme takes 
into account the long-range proximity effects. Aharonov-Bohm effects can also be included 
simply by gauging out the vector potential A separately along each lead, i.e. by multiplying 
each occurence of f m in the equations by exp(i2e A-dx/h). 

Electronic transport in a non-ideal (k,m) lead may be described by an arbitrary (but non- 
random) scattering matrix. Using the polar decomposition of this matrix, and the symmetry 
of the COE, the transmission and reflection eigenvectors can be absorbed into the scattering 
of the k and m junctions (assuming that time-reversal symmetry is preserved). This means 
that without loss of generality we may associate a transmission probability Tj with each of 
the Wk m modes in the lead. The corresponding transmission and reflection amplitudes in 



G are equal to iyTj and Jl — Tj respectively. Expressing G = G + G SG + ... in terms 
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of the TjS, fk and f m involves only inverting a j-dependent 2x2 matrix [the electron-hole 
structure is simplified by separating the even (Andreev reflecting) terms in the sum from 
the odd (non-Andreev) terms]. With the notation j3k,m = (fm — fk)/(^+fmfk)i the spectral 
current becomes 



Expressions for charge and heat transport are similarly obtained, by replacing the factor 
W k}m /\1 + fmfk\ 2 in Eqs. (|) and (|) with Y^i m {t ke ,r,me,jTtkh,j;me,j) (the minus sign is used 
for charge currents). This can be generalized to any form of Go, and so not only non-ideal 
leads but arbitrary combinations of random and non-random scattering can be treated (the 
simplification due to the polar decomposition does not easily generalize, however). 

The results are summarized in Table 1, using Nazarov's notation. gk^ m = Yjf=-T Tj 
denotes the dimensionless conductance of the lead, in the absence of proximity effects. The 
Andreev amplitudes fk are related to Nazarov's "spectral vectors" Sk by a stereographic 
projection from the complex unit disc onto the upper hemisphere: fk = tan($fc/2) exp(iipk) 
where dk and tpk are the polar and azimuthal angle of §k- The angle difference A8k, m 
is the "distance" between Sk and s m on the sphere, given by §k ■ s m = cos(A9k t m), or 
tan(A6*A;,m/2) = |/?fc, m |. In the last two lines of the table, we specialize the general expressions 
to the cases treated in Ref. || — diffusive leads and leads containing strong tunnel barriers 
(Tj — > 0) — reproducing the older results. The complex phases of the Jk, m s are also simply 
related to the directions of the spectral currents in Nazarov's picture, given by §k X s m . 

The results for (the electrical conductance of) diffusive leads can be obtained by two dif- 
ferent methods: First, by integrating over the known distribution of transmission eigenvalues 
in such leads ||12|| , one finds that the enhancement of electrical conductance by Andreev re- 
flections in relatively open channels (Tj —>■ 1) is exactly counterbalanced by the inhibition 
of "Cooper-pair tunneling" in relatively closed channels (Tj— >0), regardless of the value of 
A6k t m', Second, by treating the diffusive lead as L ^> 1 ideal (or tunnel) leads connected 
in series through L — 1 additional junctions, one finds that the spectral vectors of the in- 




(8) 
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termediate junctions lie on the big arc that connects Sk and s m with angle differences of 
A#fc m /L, resulting in very small proximity induced corrections to electrical conductivity (L 
occurences of a quantity which is quadratic in A9k, m /L). Thus, electrical conductance in 
diffusive leads is unchanged by the proximity effects being considered (see Fig. 2b below). 

Motivated in part by the "Andreev interferometer" experiment of Petrashov et al. [|13|| , 
we plot in Fig. 2 the electric and heat conductances of three structures of a geometry shown 
in the inset: one with ideal leads, one with diffusive leads, and one with tunnel barriers (each 
structure has 8 identical leads). The results are obtained numerically by repeatedly shifting 
each of the /^s in the complex plane by a fraction of the mismatch in the spectral current, 
J2 m Jk,mi until the iterations converge. An applied magnetic field leads to fast oscillations 
due to the induced phase-difference between the two superconducting electrodes, and slow 
oscillations (only one period shown) due to the Aharonov Bohm flux threaded through a 
small ring in the structure. The results exhibit some similarities to the experimental data, 
but as stated above, we are unable to reproduce the experimental oscillations (of the order 
of 10%) of the electrical resistance in a diffusive system. As pointed out in Ref. these 



may be due to the finite temperature involved (20 mK), requiring a generalization of the 
present analysis. 

In summary, using diagrammatic techniques of random scattering-matrix theory, we have 
rederived and enhanced Nazarov's circuit theory of Andreev conductance. The results follow 
in fact quite simply from notions of multiple Andreev scattering, which is taken to occur 
locally in the /cth junction with an amplitude fk (this is the amplitude of the "giant Andreev 
backscattering peak" ||). The values of fk are found self-consistent ly in a manner analogous 
to finding the voltages in the junctions of an electrical circuit. Heat transport was also briefly 
considered here, whereas the extension to finite temperatures and voltages is an important 
goal for future developments. 

I would like to thank C.W.J. Beenakker, P.W. Brouwer, B. Spivak and A. Zee for numer- 
ous discussions. Financial support of a Fulbright fellowship and NSF grants No. PHY94- 
07194 and No. DMR93-08011 is acknowledged. 

8 



REFERENCES 

[1] K. Usadel, Phys. Rev. Lett. 25, 507 (1970). 

[2] For a review, see C.W.J. Beenakker, in Mesoscopic Quantum Physics (Les Houches 
Session LXI), edited by E. Akkermans et al. (Amsterdam, Elsevier, 1995). 

[3] Yu.V. Nazarov, Phys. Rev. Lett. 73, 1420 (1994); T.H. Stoof and Yu.V. Nazarov, Phys. 
Rev. B 54, 772R (1996); See also: Yu.V. Nazarov, Phys. Rev. Lett. 73, 134 (1994); 
Report No. |cond-mat/9410011 . 



[4] N. Argaman and A. Zee, |cond-mat/9603136| , Phys. Rev. B (in press). Note that there 



the direct Andreev reflections (black semicircles) were included in G , not in S and S. 



[5] P.W. Brouwer and C.W.J. Beenakker, cond-mat/9604059 , J. Math. Phys. (in press); 
(notice esp. Sec. 6C there). 

[6] A.F. Andreev, Zh. Eksp. Theor. Fiz. 46, 1823 (1964) [Sov. Phys. JETP 19, 1228 (1964)]. 

[7] C.W.J. Beenakker, Phys. Rev. B 46, 12 841 (1992). 

[8] C.W.J. Beenakker, J. A. Melsen, P.W. Brouwer, Phys. Rev. B 51, 13883 (1995). 

[9] The conductances can be derived solely from the block off-diagonal terms in Eq. @. 
The consistency with the diagonal terms (i.e. current conservation in the leads) relies 
on the vanishing of f£ J2 m Jk,m + c.c By reversing the argument one may derive the 
expression for T, or prove that the component of Jk, m parallel to fk is conserved. 

[10] C.J. Lambert, V.C. Hui and S.J. Robinson, J. Phys. Condens. Matter 5, 4187 (1993). 

[11] The voltages V m here are the electro-chemical potentials of the quasiparticles; the elec- 
trostatic potentials are obtained by appropriately weighting this with the (vanishing) 
electro-chemical potential of the superconducting component (see Ref. |L|] below). 

[12] For a wide, diffusive wire, Y,j F(Tj) ~ g J °° dx F(cosh -2 x) for a general function F; 



9 



O.N. Dorokhov, Solid State Comm. 51, 381 (1984). Our integrations may be simplified 
considerably by a change of variables from x to 9, where tan(#/2) = tanh(x) tan(A0/2). 

[13] V.T. Petrashov, V.N. Antonov, P. Delsing, and T. Claeson, Phys. Rev. Lett. 74, 5268 
(1995); Pis'ma Zh. Eksp. Theor. Fiz. 60, 589 (1994) [JETP Lett. 60, 606 (1994)]. 

[14] Yu.V. Nazarov and T.H. Stoof Phys. Rev. Lett. 76, 823 (1996); T.H. Stoof, and Yu.V. 
Nazarov, Phys. Rev. B 53, 14496 (1996). 



10 



TABLE 



type of lead 



spectral current 



electrical conductance 



heat conductance 



ideal (Tj = 1) 

general 

tunneling (Tj — > 0) 
diffusive 



E 



W fe , m 2tan±A0 fe 

Tj sin Ae k , m 
3 1 - Tj sin 2 \Ad Km 

gk, m sinA^ 

9k,m ^@k,m 



W k , m cos 2 \A9 Km 



w COS cos # r , 

cos^ ±A6 k: „ 



cos A6^ m + 7} sin 2 \A8 Km ^ Tj cos # fc cos -d m 



^ (l-T.sin^A^) 2 
9k,m cos A9 k 

,rn 

9k,m 9k,m 



1 1 - Tj sin 2 iA0 fc , m 



g k)1n COS cos 

cos tf fc cos tf m A6V m 



sin A^ 



k,m 



Table 1: The dimensionless spectral current, | J k;m \, and electric and thermal 
conductances, I k , m /(V m — V k ) and Q k , m /(U m — U k ), for various types of leads. 



11 



FIGURES 




Fig. 1: Random-scattering-matrix diagrams for Andreev circuit theory (following Ref. [|J). 
a) Schematic drawing of a circuit with K = 3 junctions, and six external electrodes 
(M = 9), two of which are superconducting, b) Typical conductance diagram, before 
averaging; single light lines represent "deterministic" transport in the leads (Go), dangling 
double lines represent "random" scattering in the junctions, and full semicircles represent 
Andreev scattering at the superconducting electrodes, c) Couplings used to join the 
dangling lines upon averaging, and their weights, d) Diagrams for the "one-particle 
irreducible part" S, and averaged propagator G (thick line), e) Diagrams for the "two- 
particle irreducible part", T, and the "diffuson" D. 




0.5 1 0.5 1 0.5 1 



magnetic field (scaled) 

Fig. 2: Overall electrical and heat resistances between the two normal (N) electrodes of 
the circuit shown in the inset, as a function of magnetic field, in units of the resistance 
of a single lead (resistance in the absence of proximity effects = 3). a) Ideal leads; b) 
diffusive leads; c) leads with tunnel barriers. The area enclosed by the normal loop (thin 
circle) is some 30 times smaller than that of the superconducting loop (thick line in inset). 
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